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Abstract 

The modeling of carbon nanotube-metal contacts is important from both 
basic and applied view points. For many applications, it is important to design 
contacts such that the transmission is dictated by intrinsic properties of the 
nanotube rather than by details of the contact. In this paper, we calculate the 
electron transmission probability from a nanotube to a free electron metal, 
which is side-contacted. If the metal-nanotube interface is sufficiently ordered, 
we find that k-vector conservation plays an important role in determining the 
coupling, with the physics depending on the area of contact, tube diameter 
and chirality. The main results of this paper are: (i) conductance scales with 
contact length, a phenomena that has been observed in experiments and (ii) 
in the case of uniform coupling between metal and nanotube, the threshold 
value of the metal Fermi wave vector (below which coupling is insignificant) 
depends on chirality. Disorder and small phase coherence length relax the 
need for k-vector conservation, thereby making the coupling stronger. 
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I. INTRODUCTION 



Carbon nanotubes represent an intriguing new material that has attracted much atten- 
tion both from theorists and experimentalists since the early 199Os.0 Particularly exciting 
is the possibility of one dimensional metallic conductors at room temperature that can be 
used as a probe in scanning probe microscopy or as a low resistance ballistic interconnect 
for electron devices. From a more basic point of view, much can be learnt about the 
physics of conduction by studying the conductance of such a one dimensional conductor at 
low temperatures. To exploit these possibilities it is important to understand the physics of 
the nanotube-metal contacts and to experimentally demonstrate low resistance contacts in 
a reproducible manner. The contact between carbon nanotubes and metal can occur at the 
end of the tube (end-contact )@i and along the circumference of the tube (side-contact 
The low contact resistance demonstrated by de Pablo et. all and Soh et. ali are due to 
a strong interaction between metal and carbon atoms at the end of the nanotube, or/and 
due to lack of translational symmetry.! In comparison, the interaction between metal and 
carbon atoms in side-contacted nanotubes is weak. 

An interesting manifestation of weak distributed coupling is that the contact resistance 
is inversely proportional to contact length as observed experimentally by references |2] and 
|S|. Recently, Tersoff in a perceptive paper! qualitatively discussed the importance of k- 
vector conservation when the coupling between nanotube and metal is weak. The important 
physical quantities are diameter and chirality of the nanotube, Fermi wave vector of the 
metal, area of contact, and details of the metal-nanotube contact. In this paper, we study 
the physics of side-contacted nanotube-metal contactsfll by addressing how these physical 
quantities affect the transmission of electrons from the nanotube to the metal contact. For 
small diameter nanotubes, our conclusions do not fully agree with Ref. [|. We find that 
for small diameter armchair tubes, the threshold value of Fermi wave vector below which 
the conductance is very small is 27r/3ao and not 47r/3ao, which is the value for graphene. 
ao = 2.46A is the lattice vector length of graphene. In contrast to armchair tubes, the 
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threshold for zigzag tubes is zero. Our calculations also show that the conductance scales 
with contact length, a phenomena that has been observed experimentally in the work of 
Tans et. al§ and Frank et. al.@ 

In the remainder of the introduction, we discuss the salient results using simple argu- 
ments. The method is discussed in section [TT] and the numerical results and discussion are 



presented in section [TTJ. We present our conclusions in section fV . 

The first Brillouin zone of graphene touches the Fermi surface at six points (Fig. 1). Of 
these only two points are inequivalent (that is, do not differ by a reciprocal lattice vector). 
The conduction properties of graphite at low bias are controlled by the nature of eigenstates 
around these points. Consider a metal making uniform contact to graphene. The in-plane 
wave vector should be conserved when an electron tunnels from the metal to the nanotube. 
As a result, for good coupling between metal and graphene, the metal Fermi wave vector 
should be comparable to 47r/3ao, which corresponds to the Fermi wave vector of graphene. 

To discuss the case of nanotubes making contact to metal, we consider the scattering 
rate (l/r c _ m ) from the metal to nanotube within the Born approximation, 

l/r c _ m oc< ^ c \H c ^ m \ty m > , (1) 

where, ^ m (\l/ c ) is the metal (nanotube) wave function and H c _ m represents the nanotube- 
metal coupling. The wave function of an (n,m) nanotube is \l/ c ~ e tktpu (p c , where kt is the 
axial wave vector, u is the ID unit cell length, p is an integer representing the various unit 
cells and <p c is a vector representing the wave function of all atoms in a unit cell. It is 
assumed that the wave function of the metal is separable in the axial and radial directions 
of the nanotube, |\l/ m >~ e tkmVU \§ m >, where k m is the metal wave vector component along 
the nanotube axis. When the coupling between the nanotube and metal is uniform, the 
scattering rate is [Eq. flip], 

l/r c _ m oc t c . m < c |0 m > e l(km ~ kt)vu , (2) 

v 

where, the summation is performed over all unit cells making contact to metal and t c _ m 
represents a uniform coupling constant between the metal and nanotube. It is clear from 
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Eq. (|2]) that provided the metal and nanotube make contact over several unit cells, wave 
vector conservation along the axial direction is enforced as J2 P e^ km ~ kt ^ pu ~ ^5(k m — k t ). 
The axial wave vector corresponding to E — are 27r/3ao and for armchair and zigzag 
tubes respectively, and the wave vector for other chiralities varies between these two limits. 
As a result, the threshold value of Fermi wave vector below which coupling between an 
armchair (zigzag) nanotube and metal is poor is 27r/3ao (0). The threshold value of the 
metal Fermi wave vector for chiral tubes is in between that of zigzag and armchair tubes. 
As the diameter of the nanotube increases, wave vector conservation along the circumference 
becomes increasingly important, as the strip approaches a graphene sheet. 



II. METHOD 

The method used to calculate transmission probability is essentially the same as that 
in reference |], with the only addition being the connection of a metal contact .0 So in this 
section, we mainly focus on connection to the metal contact. The metal contact has a 
rectangular cross section in the (x,z) plane and is infinitely long along the y-axis as shown 
in Fig. 2. The nanotube lies on the metal contact akin to the experiment of Tans et. 
ali. In reference |2|, the nanotube bends over the edge of the metal and the influence of 
this on transport has recently been modeled by Rochefort et. al.0. In this work, the 
main focus is to model the coupling between the metal and nanotube. So we assume the 
nanotube to lie rigidly on the metal and neglect the effect of bending (Fig. 2). A perfectly 
cylindrical nanotube would touch the metal surface only along a line. To simplify modeling 
this interface, we stretch the entire circumference of the nanotube over the metal surface 
and assume coupling between carbon atoms in a sector of the circumference and the metal. 
Finally, charge self consistency^ has been neglected. 

The transmission and local density of states are calculated in a structure that can be 
conceptually divided into four parts: section of the nanotube (D), which lies on the metal 
electrode (M), semi-infinite regions of the nanotube L and R [Fig. 2]. The Hamiltonian of 



1-4 



the system can be written as, 

H = H C + H m + H c _ m and (3) 
H C = H D + H L + H R + H LD + H RD (4) 

where, H c is the pi-electron tight binding Hamiltonian of the nanotube with the on-site po- 
tential and hopping parameter between nearest neighbor carbon atoms equal to and 3.1 eV 
respectively.^ Hld and Hud are terms in the Hamiltonian coupling D to L and R respectively. 
H m and H c - m are the free particle and nanotube-metal coupling terms of the Hamiltonian. 
The Green's function G r is obtained by solving: [E — H D — Y7 L — Y7 R — E£j G r (E) = I, 
where the self energy S a = Vo a g r a V a D (a G L, R and M). g r a is the surface Green's func- 
tion of lead a and Voa {Vod) is the coupling between D (a) and a (D). The transmission 
probability between leads a and (3 \T a p] is given by, 

T aP {E) = Trace[T a (E)G r (E)T p (E)G a (E)] , (5) 

where T a (E) = 2i\Vr, a p a (E) V a n and p a (E) = —~Im[g r a (E)] is the surface density of states 
of lead a. 

The Green's function of the metal contact is calculated within the free electron approx- 
imation using the procedure outlined below. The metal contact has a rectangular cross 
section of dimensions L x and L z in the x and z directions respectively, and is infinitely 
long in the y direction. While the (yz)-coordinates are assumed to be continuous, the x- 
coordinate is assumed to be discrete with lattice spacing a = L X /(N X + 1), where N x is the 
number of lattice points. The wave functions (^? m kn) and eigen values (E m k n ) are given by, 

^mfcn(r) = X m (x)Y k (y)Z n (z) , where, (6) 
X m {x) = —=sinC^^-) , Y k (y) = -^=exp{iky) , Z n (z) = -^=sin(^^) and 

V ^X ./ Ly V L'Z 

Ti 2 rriTC h 2 k 2 h 2 titc 2 

Emkn = w^-^^TT^ + ^ + ^y ' (7) 

where, m and n are positive integers, and mo is the free electron mass. Using Eqns. (^) and 
(0) in the following equation for the Green's function, 
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g(r,r',E) = 



^ rnkny) ^ ' mkrSs ) 

E - E mkn + irj 



we obtain, 



g(r, r',E) 



im n 1 




exp[ikj\y — y'\] . mux 




h 2 L X L 



where, 



r , 9 Till o 1 r , m7T . , , l 

fc/ = i k - h- - -2 1 - cos (^rxr) + and k 




2m n E 



(9) 



For carbon nanotubes, the zero of energy (E = 0) is taken to lie at the band center. On the 
other hand, in deriving Eq. the zero of energy corresponded to the band bottom of the 
free electron metal. In the calculations, there should be only one zero of energy, which we 
take to lie at the band center of the nanotube. We also neglect charging effects and assume 
the Fermi energy of the metal to lie at the band center of the nanotube.lll Then, in the 
coordinate system where E = corresponds to the band center of the nanotube, Eq. (||) 
can be used by transforming, 



where, kf is the Fermi wave vector of the metal. 

The component of the Green's function that enters the calculation of the density of states 
and transmission probability corresponds to x = x' = a, the surface of the metal contact 
on which the nanotube lies. The (y, z) coordinates correspond to the atomic location of the 
stretched out nanotube lying on the metal. For uniform coupling between the metal and 
nanotube, we take Vdm = tD , where t is the strength of coupling between the free electron 
metal and a nanotube atom and Dq is a diagonal matrix whose dimension is equal to the 
number of carbon atoms in D. The diagonal entry Do(i, i) = 1(0) if the carbon atoms make 
(do not make) contact to the metal. 




to 




(10) 
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III. RESULTS AND DISCUSSION 



We first present results for dependence of the threshold value of the metal Fermi wave 
vector on chirality, using armchair and zigzag tubes connected to the metal contact. We 
then discuss the diameter dependence of conductance using the case of a zigzag tube as 
an example. Finally, the case of disorder in coupling between the nanotube and metal is 
considered. We consider only weak coupling between the nanotube and metal. The average 
value of the non zero diagonal elements of the coupling strength Tm are tabulated in Table 
1 for the various values of the metal Fermi wave vector considered. The main guide for the 
choice of Tm is that it be much smaller than the corresponding coupling strength between 
two carbon atoms of the nanotube (diagonal component of Tl is approximately equal to 
0.3eV for a (2,2) nanotube). A larger /smaller value of Tm results in a larger /smaller value 
of transmission in Figures 2 to 5. We calculate the transmission versus contact length 
between nanotube and metal for various Fermi wave vectors in the metal and all atoms 
around the circumference of the tube are assumed to make uniform contact with the metal. 
We emphasize that when the metal makes contact to only a sector of the nanotube such 
as in Ref. |2], the results of Fermi wave vector dependence on chirality and the conductance 
dependence on contact length are still valid. These features depend on the nanotube-metal 
coupling along the axial direction. So any change due to the finite sector will not qualitatively 
change the results. 

Experiments typically involve transmission of electrons between two metal contacts. The 
quantity Tml discussed in this section is however the transmission probability between a 
metal contact and a semi-infinte nanotube (Fig. 2). We consider this quantity because a 
long nanotube section between two metal contacts requires much more numerically intensive 
calculations. The physics discussed with regards to Tml in Figs. 3-5 hold in the case of two 
metallic contacts also, though a direct numerical comparison is not appropriate. 

In the case of armchair tubes, when the metal Fermi wave vector kf is smaller than 2ir/3ao 
(0.85A -1 ), T ML does not change significantly with contact length as shown for kf = 0.75A -1 
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in Fig. 3(a). For values of kf above the threshold, the transmission monotonously increases 
with an increase in contact length. The monotonic increase is due to weak metal-nanotube 
coupling, in which case an increase in contact length simply results in an increase in the 
transition probability to scatter from metal to nanotube.0 The transmission will eventually 
saturate with increase in contact length as there are only two conducting modes at the band 
center. For the configuration considered, Tml can have a maximum value of unity. The 
second feature of Fig. 3(a) is the increase in transmission with increase in kf. This can be 
understood by noting that electrons with a wave vector component along the nanotube axis 
that is larger than 27r/3ao scatter from the metal to nanotube, and a larger kf implies a 
large number of available metal electron states. For the purpose of these calculations, we 
considered a (2,2) armchair tube; The essential physics would in principle be true for the 
more realistic (10,10) nanotube also. 

The case of zigzag tubes is different because bands at E = cross at k = 0. Then, 
electrons in the metal electrode with any kf (no threshold) can scatter into a metallic 
zigzag tube. The results for a (3,0) tube are shown in Fig. 3(b). Here, there are two 
important points. The first point is that as there is no threshold metal Fermi wave vector, 
the transmission increases monotonically with contact length even for kf = 0.4A -1 , which is 
smaller than the threshold for armchair tubes. The second point is that the transmission for 
kf equal to 1.2A -1 is much smaller than that for armchair tubes [Fig. 3(a); the transmission 
of the three smaller values of kf have been multiplied by a factor of ten.]. This is because the 
nanotube wave vector around the circumference (k c ) of a zigzag tube is large; k c = 4ir/3ao for 
the crossing bands and as a result, the overlap integral [Eq. ([!])] is small. As kf = 1.75A -1 is 
larger than the threshold for graphite, the transmission probability is larger, and comparable 
to that for armchair tubes [Fig. 3(b)]. 

What happens when the diameter increases? In the limit of large diameter, a nanotube 
is akin to graphene and the threshold kf to couple well with metal should approach 47r/3a J 
Numerically, it is difficult to simulate large diameter tubes along with large contact lengths 
because of time and memory requirements associated with the calculation of g r M . To convey 
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the main point we consider two simpler cases, the first case compares the transmission 
probability of the two smallest semi-metallic zigzag tubes with varying contact lengths and 
the second case considers zigzag tubes of varying diameters with a rather small contact 
length. Fig. 4 compares the transmission probability versus contact length of the (3,0) 
and (6,0) nanotubes; The (6,0) nanotube has double the diameter of the (3,0) nanotube. 
The (6,0) nanotube correspondingly has a smaller transmission and the trend of decrease in 
transmission will continue with further increase in diameter. The inset is a calculation of 
transmission probability versus diameter of semi-metallic zigzag tubes for a contact length 
of 42.6 A (ten unit cells). Tml decreases with increase in diameter because wave vector 
conservation becomes increasingly important with increase in diameter. Shown also in this 
figure for comparison are 1 / diameter and 1 / \J diameter. 

We now address the role of disorder. Disorder in either the nanotube, metal or nanotube- 
metal coupling will in general result in larger transmission when compared to the disorder- 
free case. Wave vector conservation is relaxed due to scattering from defects and transmission 
will increase with increase in contact length even when the metal kf is below the threshold 
value. We consider the case of disorder in nanotube- metal coupling (H c _ m ). Disorder in all 
elements of the coupling between the nanotube and metal was introduced randomly. The 
disorder in coupling of atom i to the metal contact can be written as, U = at av + (l — a)t™ nd , 
where t av is the average value of U over all sites connected to the metal and a is a fraction 
between zero and unity. t\ and is the random component whose average is equal to t av . In Fig. 
5, the two strengths of disorder correspond to a = and a = 0.5 (smaller a corresponds 
to larger disorder), such that t av has the same value as in Fig. 3(a). For an armchair tube 
in contact with a metal with kf = O.TSA^ 1 , the transmission was very small and more 
importantly did not vary with contact length [Fig. 3(a)]. Introducing disorder changes 
this trend and causes a monotonic increase in transmission with length of contact [Fig. 5]. 
Similarly, for large diameter tubes, in the presence of disorder there should be significant 
transmission when kf is smaller than the threshold 47r/3ao- The requirement of wave vector 
conservation is also relaxed when the phase coherence length is small. So we expect the 

1-9 



coupling to improve with decrease in phase coherence length. 

IV. CONCLUSIONS 

In this paper, we addressed some aspects of the physics of a nanotube side-contacted to 
metal, a problem of current importance. Coupling of carbon nanotubes to metal depends 
on both chirality and diameter. Wave vector conservation of an electron scattered from 
the nanotube to metal plays a central role in determining the transport properties. The 
difference between small and large diameter nanotubes is that while in the former wave 
vector conservation is important only in the axial direction, in the latter it is important 
in both the axial and circumferential directions. As a result, small diameter armchair and 
zigzag tubes have a cut-off value of the metal Fermi wave vector equal to 27r/3ao and zero, 
respectively. For chiral tubes, the cut-off value of the metal Fermi wave vector lies in between 
these two limits, with the value decreasing with increase in chiral angle. A large diameter 
nanotube is akin to a graphene sheet and the cut-off value of the metal Fermi wave vector 
in this case approaches 47r/3ao with increase in diameter. Disorder in the metal, nanotube 
or their coupling relaxes the requirement of k-vector conservation and in general improves 
coupling. The groups of references |] and || have shown increase in conductance with contact 
length. In this paper, we discussed two situations that could lead to this. The first situation 
requires the metal Fermi wave vector to be larger than the threshold discussed in the text 
and holds even when there is no disorder. The second situation requires disorder in coupling 
to the metal but there is no restriction on the value of the Fermi wave vector. 
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Figure Captions: 

Fig. 1: First Brillouin zone of graphene. Points P, P', P", Q, Q', Q" touch the Fermi 
surface. a is the lattice vector length of graphene. A metal with Fermi wave vector smaller 
(inner circle) and larger (outer circle) than 47r/3ao couples poorly and well to graphene 
respectively. AE NC is the energy difference between the first non crossing subband below 
and above E=0. 

Fig. 2: A metal making contact to a nanotube. The (x, z) dimensions of the metal form a 
rectangular cross section with lengths (L X ,L Z ). The y direction is infinitely long. 

Fig. 3: Transmission probability for (a) armchair and (b) zigzag tubes versus contact 
length. In both cases the largest contact length corresponds to sixty unit cells. The main 
point of (a) is that for the metal Fermi wave vector smaller than the threshold 27r/3a , 
coupling between the nanotube and metal is small and increasing the contact length does 
not change the transmission probability. For metal Fermi wave vector larger than 2n/3ao, 
the transmission probability increases with increase in contact length and also with increase 
in kf for a given contact length. The main point of (b) is that there is no threshold in the 
metal Fermi wave vector. Even in the case of a small value of the metal Fermi wave vector 
(0.4A -1 ), the transmission increases with increase in the contact length, albeit the magnitude 
of transmission is small. As in the armchair case, the transmission probability increases with 
increase in kf for a given contact length. The values of T ML in (b) corresponding to kf equal 
to 0.4, 0.75 and 1.2A are multiplied by a factor of ten. 

Fig. 4: Comparison of transmission probability of (3,0) and (6,0) nanotubes versus contact 
length. The transmission probability decreases with increase in diameter. Inset: The y-axis 
is T M l for metallic zigzag tubes scaled by 1.0 e+4. The solid line is the diameter dependence 
of T ML for a contact length of 42.6 A. The upper and lower dashed lines are 1/Vdiameter 
and 1/diameter dependences, shown for comparison. 
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Fig. 5: Comparison of transmission probability versus contact length for a (2,2) armchair 
tube, with and without disorder in nanotube-metal coupling. The metal Fermi wave vector is 
0.75A -1 . Note that for the case without disorder, the transmission is poor and increasing the 
contact length does not help. Introducing disorder changes this picture and the transmission 
begins to increase with increase in contact length because k-vector conservation is relaxed. 
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TABLES 



Table 1: 



Metal Fermi 
Wave Vector 

(A" 1 ) 


r M (e V) 


0.4 


1.1 x 10" 4 


0.75 


7.2 x 10" 4 


0.9 


1.2 x 10" 3 


1.2 


2.9 x 10" 3 


1.75 


9.1 x 10" 3 



Tm for the different values of the metal Fermi wave vectors used. 
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FIGURES 
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Fig. 2 / Anantram 
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Fig. 3 / Anantram 
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Fig. 4 / Anantram 
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Fig. 5 / Anantram 
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